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^D i Abstract. According to previous investigations the effect of diffusion in the 

^N ' stehar atmospheres and envelopes of subdwarf B (sdB) stars with luminosities 

10 <^ L/ Lq <^ 100 strongly depends on the presence of weak winds with mass 

loss rates M <^ IO^^^Mq/yi. These calculations with the mass loss rate as a 
free parameter have shown that it is hardly possible to reproduce the measured 
abundances of helium and metals simultaneously. A possible reason is the decou- 
pling of metals, which preferably absorb the photon momentum, from hydrogen 
i-G ■ and helium in the wind region. In the present paper it will be investigated if 

Mh' "chemically homogeneous" winds, as assumed in previous investigations, with 

Q . mass loss rates M < 10~^^MQ/yr can exist. From an observational point of 

Jh ' view the existence of weak winds in sdB stars is unclear. Only in the most 

c/j I luminous ones possible wind signatures have been detected. Therefore it will be 

C^ . investigated if according to the theory of radiatively driven winds the existence 

of weak winds is plausible. A stellar mass Af» = Q.5Mq is assumed. 

The resuhs for effective temperatures Teg = 35000, 30000 and 25000K, 
^ ■ metallicities 0.1 < Z/Zq < 1 predict decreasing mass loss rates with increasing 

(^ ' surface gravity. Dependent on the luminosity and metallicity the mass loss rates 

Q^ , are between about IO^^^Mq/yi and zero. If at all, chemically homogeneous 

^^ ■ winds can exist for the most luminous sdB stars only. For the other ones selec- 

tive winds are expected which should lead to additional changes of the surface 
r — ^ , composition. 

(^ ■ In sdB stars, hot white dwarfs and HgMn stars (which are chemically pecu- 

(yQ I liar main sequence stars) the measured metal abundances are tendencially lower 

f^ , than the ones predicted from diffusion calculations which assume an equilibrium 

between gravitational settling and radiative levitation. Only for helium in al- 
most all cases the measured abundances are larger than the predicted ones, but 
usually lower, below the solar value. This may be an indication that the abun- 
^ ■ dance anomalies of metals are preferably due to the selective winds, whereas 

Cu I the helium deficiencies are due to gravitational settling, which for still unknown 

reasons is less effective than expected in an undisturbed stellar atmosphere. 



1. Introduction 

The abundance anomalies in subdwarf B (sdB) stars, white dwarfs and chemi- 
cally peculiar main sequence stars are believed to be at least partially due to the 
effect of diffusion in the stellar atmosphere and envelope. Several attempts have 
been made to explain the abundances with the effect of gravitational settling 
which may be counteracted by radiative levitation. As the radiative force on 
an element decreases with increasing abundance due to saturation effects, the 
surface composition can be predicted from the equilibrium condition between 

1 






Unglaub 



the inward gravitational force and the outward radiative force. However, an 
agreement between predicted and measured abundances requires the absence of 
disturbing processes Uke mass loss or convective mixing. This may be the rea- 
son why in many cases the agreement is not satisfactory. A com parison between 
jredicted and measured metal abundances for sdB stars f e.g. [Bergeron et al 



predicted and measured metal abundances tor sdo stars l e.g. [.Bergeron et ai. 
[l988 : Ohl et al.l l2000l : IChaver et al. 2006: iBehara fc Jeffer"vl Ithese DroceedingsT) 



and hot white dwarfs (e.s. ISchuh et alTl200l l2005l : iGood et al.1 l2005l ^ shows 
that the measured abundances are tendencially lower than the predicted ones, 

although a few exceptions exist (e.g. silicon in hot white d warfs). 

According to recen t spectral analyse s of s dB stars (e.g. Geier et al.ll2008al lbl: 
O'Toole k Heberll2006l : lEdelmann eralll2006l : iBlanchette et al.ll2006l ) esoeciallv 



the hotter ones with Tgff > 30000K in many cases show strong deficiencies of 
light metals like Al, Mg, O and Si by more than a factor of 100 in comparison 
to the solar abundances, whereas enrichments of elements heavier than iron 
by a factor of 100 are not unusual. Qualitatively, these abundance patterns 
show some similarity to those found in HgMn stars, whic h are a subg roup of 
the chemically peculiar main sequence stars reviewed by ISmithI ( 19961 ). The 
HgMn stars with lOOOOK ^ T^s < 16000K, log 5 ^ 4.0 are characterized by 
low rotational velocities and weak or non-detectable magnetic fields. Some light 
metals (e.g. Al, N) tend to be deficient, whereas heavy metals (e.g. Hg, Mn, Pt, 
Sr, Ga) may be enriched up to several orders of ma gnitude. So me of the mos t 
recent spectral analyses are from lZavala et ahl ( 2007) an d Adelma n et al. ( 20061). 
The measured abundances of iron group elem ents ( Seatonlll996l : IJomaron et al.l 
19991 ) as well as nitrogen ( Robv et al.lll999l ) tend to be lower than predicted 
from equilibrium diffusion calculat ions. Only for mer cury abundances larger 



than predicted have been detected ( Proffitt et al.|[l999l ). 

So in hot white dwarfs, sdB and HgMn stars a common tendency seems to 
be present, according to which the metal abundances are lower than expected 
from the equilibrium condition between gravitational settling and radiative levi- 
tation. In addition there is a large scatter of abundances from star to star. Even 
for stars with similar stellar parameters the abundances may be different. This 
points to some time-dependent process and not to an equilibrium state. In main 
sequence stars the chemically peculiar phenomenon is restricted to stars with low 
rotational velocities. In addition the presence of magnetic fields may be of impor- 
tance, beca use magnetic fie lds may change the radiative acceleration or suppress 
convection ( Turcottg l2003l ) . White dwarfs and sdB stars, however, are always 
more or less chemically peculiar. Up to now no correla tion between abund ance 
anomalies and magnetic field strengths has been found ( Q'Toole et al.ll2005l ). In 
hot DA white dwarfs no outer convection zones should exist, because hydrogen 
is preferably ionized and helium is strongly deficient. In sdB stars a thin super- 
ficial convection zone with a mass depth of the order W~^'^M* may be pr esent 
only for helium abundances He/H > 0.01 by number ( Groth et al.lll985l ). So 
it seems to be unlikely that magnetic fields are of decisive importance for the 
explanation of the surface compositions. 

In hot DAO white dwarfs helium is detec table, but in man y cases it is de- 
ficient in comparison to the solar value (e.g. Nap iwotzkJ^19 99|). The sam e is 
true for the majority of sdB's (e.g. lEdelmann et al . 2003; Liske r et al.ll2005l ) and 
the helium deficient main sequence stars (e.g. in HgMn stars helium usually is 
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deficient). In contrast to the metal abundances, however, the abundances of 
hehum are alw ays larger than predicted from equilibrium diffusion calculations 
( Vennes et al.i 1988 : Michaud et aUll989l , [l979 ). These results may be explained 
with the effect of gravitational settling which, however, somehow must be dis- 
turbed. Possible reasoi is for this disturba nce of the equil ibrium may be the 
presence of turbulence (Vauclair et al.lll978i ) or mass loss. iFontaine Sz Chavei 
( 19971 ) and lUnglaub &: BuesI ( 19981 ) predicted helium abundances as a function 
of time in the presence of weak winds. The results have shown that for mass loss 
rates of the order lO~^^M0/yr helium sinks much more slowly than in the case 
of an undisturbed stellar atmosphere. Within the lifetimes of sdB stars near the 
extended horizontal branch (~ lO^yr) the helium abundance would gradually 
decrease from the solar value to He/H « 10~^ by number. This could explain 
why the helium abundances usually are in this range. 

If this scenario with weak winds were the correct explanation for both the 
helium and the metal abundance anomalies, then it should be possible to find a 
mass loss rate for which all abundanc es can be explained simultaneously. The 
calculations of lUnglaub Sz Bued ( 200ll ) for the elements H, He, C, N and O have 
shown that this is hardly possible. According to these calculations for solar 
initial composition helium should always be more deficient than the metals. No 
mass loss rate exists which leads to deficiencies of C and O by more than a factor 
of 100, whereas helium is deficient by a fac tor of ten only. Th is, however, is not 
an unusual composition in sdB stars (e.g. iHeber et al.ll2000l ). Moreover in the 
presence of winds with mass loss rates of the order lO~^^M0/yr which are re- 
quired to explai n the helium abu ndances, the proposed pu lsation mechanism of 
some sdB stars ( Charpinet et al.|[l9 97: Fo ntaine et al.ll20C)3l ) would become ques- 
tionable. As mass loss tends to level out concentration gradients, the reservoir 
of iron (or other iron group elements) in the stellar envelope needed to explain 
the pulsations should be destroyed in t ime scales wh i ch are mu ch shorter than 
the lifetimes of sdB stars. According to lChaver et al.l (J2004l ) and lFontaine et ah 
( 20061 ) for a mass loss rate of M = 6 x lO~^^M0/yr the reservoir would be 
destroyed on about lO^yr. For M = lO~^^M0/yr the matter in mass depths 
^ 10~^M*, where the reservoir is expected, would be blown away within one 
million years only. 

Probably the most crucial assumption in these diffusion calculations with 
mass loss has been that the winds are "chemically homogeneous" . If Mi is the 
mass loss rate of an element I and Q_ its mass fraction in the photosphere, then 
this assumption states that Mi = QM, where M is the total mass loss rate. Such 
a chemically homogeneous wind prevents (if M ^ lO~"'^"^M0/yr) or retards (if 

M < lO^"'^^M0/yr) gravitational settling. However, it does not directly change 
the surface composition. The opposite case would be a "selective" wind in 
which the mass loss rates of the individual elements are essentially independent 
of each other. A selective wind should lead to additional changes of the surface 
composition, which have not yet been taken into account in the calculations. 

In radiatively driven winds of hot stars the photon momentu m is absorbed 
preferably by the metals (see e.g. lAbbottlll982l : IVink et aUl200ll ). whereas the 
contribution of hydrogen and helium is small. Thus the metals are acceler- 
ated and move outwards. If the flow of metals is sufficiently large, then due to 
Coulomb collisions with the metals hydrogen and helium are accelerated as well. 
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For this purpose, in dense winds a small velocity difference between the metals 
and hydrogen and helium is sufficient. Then it should be a good approximation 
that all elements have the same velocity and that the wind is chemically homo- 
geneous. If, however, the flow of metals is small, then it may happen that the 
coupling of the various constituents due to collisions is not sufficiently effective 
and that hydrogen and helium are left behind. Th is scenario may lead to pure 
metallic winds such as investigated by iBabell ( 19951 ) for main sequence A stars. 

Mass loss has b een detected in subdwarf stars (sdO) which are more 
luminous than sdB's ( Hamann et al.lll98ll : lRauchlll993l ). Up to now there is no 
observational proof for the existence of wind s in sdB stars. From a quantitative 
analysis of Ha line profiles of 40 sdB stars ( Maxted et al.ll200ll ). a comparison 
of synthetic NLTE Ha line profiles from static model atmospheres with the 
observations revealed perfect matches for almost all stars. Only in the four most 
luminous sdB's anomalous Ha lines with a small emission at the line center hav e 
been detected, which possibly are signatures of weak winds (JHeber et al.ll2003l ). 
For the case T^q = 36000K, log 5 = 5.5 and l og L/Lfn = 1.51, from a spectral 
synthesis of Ha with his wind code lVinkI ( 2004i ) found a similar behaviour of Ha 
if the existence of a weak wind with M ~ 10~^"'^MQ/yr is assumed. 

In Sect. 2 it will be investigated for a stellar mass M* = 0.5Mq if winds 
with mass loss rates M <^ IO^^'^Mq/jt can be chemically homogeneous. The 
argument s are similar as i n the investigations for more l umin ous stars e.g. 
fromOwocki Sz Pulsl ( 20021 ) . ISpringmann fc PauldrachI (j 19921 ) and lKrticka et al. 
( 20031 ). All metals are lumped together into one mean metal which is acceler 



at ed due to the absorption of photon momentum. This is some simplification, 
in iKrtickal ( 20061 ) the individual elements are considered separately. Hydrogen, 
helium and the free electrons are denoted as "passive plasma" which can only be 
accele rated due to co llisio ns with the outflo wing metals. From the calculations 
e.g. of I Abbott] ( 19821 ) and lVink et al.l ( 200ll ) as well as from own calculations as 
described in Sect. 3 there is no indication that the radiative force on hydrogen 
and helium is not negligible. The mean radiative acceleration on the metals 
exceeds the one on hydrogen and helium by at least a factor of 100. 

In Sect. 3 the results of mass loss calculations for sdB stars according 
to the origin al theory of radiatively driven winds of ICastor. Abbott. &: Klein 
(|l975. CAKI ) are p resented and ar e corn pared with the predictions from the 
mass loss recipe of IVink &: Cassisil ( 20021 ). For these wind models, which are 
obtained from a one component description of the wind, it is again checked if 
the metals may be coupled to hydrogen and helium. In Sect. 4 the consequences 
of the results for the surface composition of sdB stars are discussed. 



2. Do Chemically Homogeneous Weak Winds Exist? 



In the following the hydrogen and helium, for which the radiative acceleration 
is assumed to be zero, will be denoted as "element" 1, whereas the mean metal 
is denoted as "element" 2. As in the supersonic region the gradient of the gas 
pressure can be neglected, the momentum equation for "element" 1 (H+He) can 
be written as: 

ffcou = — 9- + vi-r- (1) 

r^ ar 
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gcoU is the collisional acceleration on hydrogen and hehum due to Coulomb 
interactions with the metals and vi is the mean velocity of hydrogen and helium. 
The existence of a supersonic wind with increasing velocity in outward direction 
requires that the acceleration term vi-^ is larger than zero. Thus this equation 
can only be valid if: 

ffcoU > — — (2) 

where G is the gravitational constant and r the radius. Hydr ogen and helium 



wnere G is tne gravitational constant and r the radius, nydr ogen and helium 
cannot be accelerated if this condition is violated. As derived by lBurgersI ( 19691 ) . 
^fcou can be written as: 

ffcoii = r;^; (In A) G(^) (3) 

mim2 kl ^ ' 

rrii and m2 are the mean masses of "elements" 1 and 2, p2 is the mass density 
of the mean metal, k the Boltzmann constant, T the temperature and e is the 
electron charge (cgs system). Zi is the mean charge of hydrogen and helium, Z2 
is the charge of th e mean metal. In A is the Coulomb logarithm which according 
to lBurgera (|l969l ) is: 



lnA = -l/2 + lnft^V (4) 



\Z1Z2ey 
Rb is the Debye radius: 



0.5 



/ kT 

^° " (,47re2(niZ2 + n2Z| + ne)j ^^^ 

where ni is the number density of hydrogen and helium, n2 and Ug are the 

number density of the mean metal and the free electrons, respect ively. 

Gfx) in Eq. 3 is the Chandrasekhar function defined as, e.g., in lKrticka et al. 



( 20031 ). It depends on the quantity x, which approximately is the velocity differ- 



ence of the mean metal and the passive plasma in units of the thermal velocity 
of hydrogen. If V2 is the velocity of the mean metal and vi the mean velocity 
of the passive plasma, then the latter can be accelerated only if ^2 — ^i > 0. 
For small velocity differences, Gm and thus the collisional acceleration approx- 
imately is proportional to t'2 — tii. If, however, the velocity difference is near 
the thermal velocity of hydrogen, Gm reaches a maximum value and decreases 
to larger velocity differences. For the discussion in the present section only this 
maximum value is important: 

G^ax = 0.214. (6) 

From Eq. 3 it can be seen that the collisional acceleration on hydrogen and 
helium is proportional to T~^. In the present paper T = Tes is assumed for 
the wind region. If the temperature were larger, e.g. due to frictional heating, 
this would reduce ^coU and thus favour the decoupling of the metals. This effect 
could (partially) be compensated, if the heating leads to a higher mean charge 
of the metals, because (7coii ~ -^1 ■ Thus in this case more detailed calculations 
would be necessary. 



6 Unglaub 

The mass density p2 of the mean metal in Eq. 3 can be substituted from 
the equation of continuity for the metah 

M2 ^^^ 

47rr^f2 

where M2 is the mass loss rate of the metals and V2 their mean velocity. With 
this substitution condition Eq. 2 can be transformed into a condition for M2: 

This condition depends on the mean velocity of the metals. It is most restrictive 
in the part of the wind where f 2 is maximal. In accelerating winds this maximal 
velocity is the terminal velocity foo- Ac cording to theoretical calc ulations and 
observations of hot star winds (see e.g. iLamers &: Cassinellilll999l ). radiatively 



driven winds accelerate to terminal velocities which are at least of the order of 

the surface escape velocity Ugsc- Thus we require that this condition must still 

be fulfilled for 

/2GM* 1 

V2 = Vesc = J -^ = V2 {GM,g) 3 . (9) 

Here it has been assumed that the radiative force due to electron scattering is 
negligible which is justified for subluminous stars. The right equation is obtained 
from Rf, = y/GM^jf^. With this expression for V2 condition (8) can be written 
as: 

Now we assume a mean mass for hydrogen and helium mi = nip (where nip is 
the proton mass) and ni2 = 15nip for the mean metal. The mean charges are 
assumed to he Zi = 1 and Z2 = 3, respectively. According to own calculations 
In A may vary between about 6.0 at the wi nd base and 18. in th e outermost 
regions. As an upper limit and similar to lOwocki fc Puld ( 20021 ) we assume 



InA = 20. With a typical stellar mass for sdB stars M^, = 0.5Mq and for 
T = Teff it follows: 

M2 > 1.2 X lO^^^Teff^i . (11) 

The mass loss rate M2 of the metals is given in Mq/ji, Tgff in K and the 
surface gravity g in cms~^. For typical stellar parameters of sdB stars, e.g. 
Teff = 30000K, logg = 5.7, it follows 

Ma^lO-i^Mo/yr. (12) 

This means that the mass loss rate of the metals alone must exceed this value. 
Otherwise hydrogen and helium will decelerate and cannot be expelled from the 
star, if not the decoupling occurs at such a large radius at which the velocity is 
already larger than the local escape velocity. 

If the wind is chemically homogeneous, this implies M2 = C2M. If the 
mass fraction of the metals is of similar order of magnitude as the solar one 
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(C2 ~ 0.01), then the condition M2 ^ lO~^^M0/yr implies for the total mass 

loss rate M ^ lO~"'^'^M0/yr. Thus chemically homogeneous winds with lower 
mass loss rates cannot exist, if the met allicity is not larger tha n solar. A similar 
result can be obtained from eq. 22 of lOwocki Sz Pula ( 20021 ). In this paper a 



maximum velocity Umax is derived up to which the constituents may be coupled. 
If the wind accelerates at least to the surface escape velocity, then a necessary 
condition for a the existence of a coupled wind is Vmux > Ve 



^esc 



3. Mass Loss Rates from a One Component Description of the Wind 

In this section it will be discussed how large the mass loss rates are in sdB stars 
according to the theory of radiatively driven wind s. These calculations will b e 



described in more detail in a forthcoming paper (jUnglaubl lAfcA. submittedl ). 
The momentum equation for an isothermal wind is solved without the usual 
parametrization of the l ine force multip lier parameters, which becomes ques- 



tionable in weak winds ( Kudritzkil I2OO2I ) . As later improvements of the CAK 



theory (finite disk correction and changes of ionization in the wind) are ne- 
glected, the method of solution is straightforward. The omission of the finite 
disk c orrection leads to an ove restimate of M by a factor of the order two to 



three (Lamers &: CassinellilllQQS). From these calculations, in which multicom- 



ponent effect are neglected, the total mass loss rate M and the velocity law U(^) 

are derived. From this solution the acceleration term u^ is known. With the 
assumption that metals are trace elements {vi ~ v) it can be checked if met- 
als may be coupled to hydrogen a nd helium from a criterion similar to the one 
derived by lOwocki &: Pula ( 20021 ). This criterion requires that hydrogen and 



helium must accelerate. Thus it will predict decoupling at larger mass loss rates 
than the less restrictive criterion Eq. 2, which only requires that the passive 
plasma does not decelerate. 

In Fig. 1 the predicted mass loss rates are shown as a function of surface 
gravity for Tefj = 35000, 30000 and 25000 K and several metallicities: Z/Zq = 1, 
1/ 3 and 1/10. For ZlZ p) = 1 and 1/10 the results are compared to the ones 
of lVink fc Cassisil ( 20021 ). which are represented by dotted lines. It can be seen 



that especially for Z/Zq = 1 the agreement is very well in many cases. This, 
however, is a consequence of the assumptions. As in the present calculations the 
original version of the CAK theory is used, for the terminal velocity it follows 
v^ ~ Uesc and the solution of the momentum equation has the form 



T R* 

'J(r) =V^yl- — 



0.5 



Vink Sz Cassisil ( 20021 ) assumed a similar velocity law with foo as a free param- 



eter. Their results shown in Fig. 1 are for Uoo = Wesc- With this v^r) they obtain 
the mass loss rate from the requirement of a global momentum conservation. As 
the function Vfj.\ approximately is the same in both calculations, the agreement 
of the mass loss rates should be expected. Nevertheless it may appear surpris- 
ing, because Vink & Cassisi's calculation of the radiative acceleration is clearly 
more sophisticated than the present one. They took into account about 10^ lines 
of the elements H~Zn with NLTE occupation numbers for the most important 
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Figure 1. Predicted mass loss rates as a function of surface gravity for 
Toff = 35000, 30000, 25000K and Z/Zq = 1, 1/3, 0.1 with M* = 0.5Mq. 
Dashed Unes indicate decouphng of metals in the corresponding wind model. 
The upper and lower dotted lin es in each figure represent the results of Vink 
& Cassisi (iVink fc Cassisil[200l) . 
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elements, whereas in the present calculations only the elements H, He, C, N and 
O have been taken int o account with a line list similar to the one used in the 
diffusion calculations of lVauclair et alj ( 19791 ) . It consists of about 150 lines only 



and the occupation numbers have been obtained with the assumption of LTE, 
which is not a good approximation at low densities in the wind region. 

The main reason why in spite of these very different assumptions the pre- 
dicted mass loss rates are similar is, that in we ak winds the major part of th e 
radiative force is due to a few strong lines only (JAbbottlll982l : |Puls et al.ll200Cll ). 



E.g. for Tes = 30000K, logg = 5.5, Z/Zq = 1 the predicted mass loss rate ac- 
cording to the present calculations is M = 4.7 x lO^^^M^g/yr. In this example, 
about 70 % of the radiative acceleration are due to three carbon lines only: these 
are the hues C III A977 A and Al 176 A and C IV A 1549 A (doublets or multiplets 
are considered as one line). If only these three lines were taken into account, the 
mass loss rate would be lower by a factor of two only. This means that the result 
essentially depends on the abundance of carbon. With the present assumptions 
58 % of ah carbon is C III and 42 % is C IV. If the C III/C IV ionization equi- 
librium is changed, then the mass loss rates would still be of similar order of 
magnitude. With the assumption that all carbon is in the ground state of C III 
so that from the considered three lines only A977A contributes, the predicted 
mass loss rate would be lower by a factor of two. If all carbon were in the ground 
state of CIV, so that only A1549A contributes, M would decrease by a factor 
of three. So independently how this ionization equilibrium is, at least the order 
of magnitude of the mass loss rates will always be similar. 

The fact that the wind is driven by a few strong lines may have impor- 
tant consequences, however. In the CAK theory it is assumed that the stellar 
flux which can be absorbed by the matter in the wind is independent of the 
velocity. This is justified in dense winds which are driven by a large number 
preferably of weak l ines. In the present calc ulations the emergent flux of models 
as described e.g. in lUnglaub fc Buea ( 200ll ) has been used with a temperature 



structure according to the diffusion approximation for the flux. This flux is in 
good agreement with the continuum fluxes of model atmospheres of sdB stars 
(Heber, priv. comm.). Near the center of strong lines, however, the flux may be 
lower by a factor of the order ten. Thus in the inner parts of the wind, where 
the velocity and thus the Doppler shift is small, the present calcuations clearly 
overes t imate the stellar flux and thus the radiative acceleration. According to 
Babell ( 19961 ) this effect of "line shadowing" leads to lower mass loss rates and 



larger terminal velocities. So the present calculations probably overestimate the 
mass loss rates. This may also be true for Vink & Cassisi's results shown in 
Fig. 1. With the assumption Uoo = 5uesc instead of Voo = Vesc, their mass loss 
rates would be lower by a factor of the order of ten. 

The predicted mass loss rates shown in Fig. 1 clearly depend on the surface 
gravity and the metallicity. For Z/Zq = 1 coupled winds can exist only if log g ^ 
5.5. As the criterion for ion decoupling is more restrictive than the one used in 
Sect. 2, this requires total mass loss rates almost of the order lO~^^M0/yr. If 
the metallicity is reduced by a factor of ten, then for the considered effective 
temperatures at least for log g > 5.5 no solution of the momentum equation for a 
one component plasma exists at all, because the radiative force is not sufficient. 
This result has been confirmed by calculations of Vink (priv. comm.). For sdB 
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Figure 2. Lines in the Teff — log 5 diagram above which chemically homoge- 
neous winds may exist for Z/Zq = 0.1, 1/3 an d 1, respective l y. Sq uares 
and circles r eprese nt the sdB stars analyzed by iMaxted et al.1 ()2001[ ) and 
iLisker et al.l (|2005| ). respectively. Filled symbols represent the sdB's with 
peculiar Ha line profiles, which may indicate the presence of a weak wind. 



stars with log 51 = 5.5 and Z/Zq = 1/10 he could not find any mass loss rate for 
which a global momentum balance is fulfilled. 



4. Discussion 

In the Teff - log 5 diagram of Fig. 2 for 25000K < Tcs < 40000K and for var- 
ious metallicities the lines are shown above which according to the mass loss 
calculations as described in Sect. 3 chemically homogeneous winds may exist. It 
can be seen that for the majority of sdB stars this is not possible if the metal- 
licity is solar or subsolar. They are below the line for Z/Zfr, = 1. T h e sdB 
star s introduced in th e diagram have been analyzed by iMaxted et al.l ( 200ll ) 
and iLisker et al.l ( 20051 ). According to the assumptions of the present paper in 
those ones with peculiar Ha line profiles (represented by filled symbols) chemi- 
cally homogeneous winds may indeed exist if the metal abundances are not too 
far below the solar value. The predicted mass l oss rates agree with the ones 
from the mass loss recipe of lVink &: Cassisil ( 20021 ) and are of the order 10"^" to 
lO~^^M0/yr. However, as explained in Sect. 3, this agreement does not exclude 
that the mass loss rates are overestimated. So this result is still questionable. 
The existence of winds may depend on the abundance of one element only, which 
preferably contributes to the radiative acceleration. 

From the results it is clear that chemically homogeneous winds with mass 
loss rate s M <; lO~^^-M0/yr c annot exists. From arguments similar as in the 
paper of lOwocki &: Pulsl ( 20021 ). this result can be obtained without the calcula- 
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tion of a wind model. The mass loss rate of the metals alone must be at least of 
the order lO~^^M0/yr. Otherwise hydrogen and helium cannot be accelerated 
throughout the wind. This result could be questioned only if a wind solution 
could be found for which the terminal velocity of the metals is clearly lower than 
the surface escape velocity. iKrticka &: Kubatl (J200 (j) suggested that the termi- 
nal velocity could indeed be lower than Uesc (by a factor of the order two only, 
however) , if the wind switches to a "shallow" solu t ion w ith a n abrupt lowering 
of the velocity gradient. However, lOwocki &: Pulsl ( 20021 ) and lKrticka Sz Kubatl 



( 2OO2I ) argued that these solutions are unstable. 



If hydrogen and helium cannot be expelled from the star, then pure metallic 
winds may exist. As the outflowing metals in the stellar atmosphere not only 
have to overcome the gravitational force, but in addition the frictional force due 
to collisions with protons and helium particles, the metal abundances should 
be lower than predicted from equilibrium diffusion calculations (if concentra- 
tion gradients are negligible). For these metals, for which the mass loss rate is 
sufficiently small, measured and predicted abundances should be in agreement. 
For metals with non-zero mass loss rate this scenario shou l d lead to abundances 
varying with time as has been discussed bv ISeatonI ( 19961 . Il999l ) for iron group 
elements in the envelopes of HgMn stars. 

If both hydrogen and helium are in hydrostatic equilibrium, then measured 
helium abundances should be approximately in agreement with the ones pre- 
dicted from equilibrium diffusion calculations. The fact that in the various types 
of helium deficient chemically peculiar stars the measured ones are larger, seems 
to indicate that gravitational settling in general is less effective than expected 
in an undisturbed stellar atmosphere. Moreover the existence of two distinct 
sequen ces of sdB stars which are characterized b y an offset in the helium abun- 
dance ( Edelmann et al]|2003l : iLisker et al.ll2005l ) can hardly be explained with 



one atmospheric effect alone. It may point to a dependence on the star's history. 
Several scenarios of single star and binary evolutio n of the sdB and the hotter 
sdO stars are under discussion (jStroeer et al.l 120071 ). 
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